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ABSTRACT 

Recent observational studies have revealed star-to-star abundance inhomogeneity 
among light elements (e.g., C, N, O, Na, and Al) of stars on the main sequence in the 
Galactic globular clusters (GCs). One of promising interpretations for this result is 
that the observed abundance inhomogeneity is due to the second generation of stars 
formed from ejecta of the first generation of evolved stars (e.g., AGB stars) within GGs. 
However it remains unclear whether and how this primordial pollution can occur within 
GCs. We here propose a new scenario in which primordial pollution of GCs is highly 
likely to occur if GCs are located in the central regions of high redshift dark matter 
subhalos that can host low-mass dwarf galaxies. In this scenario, gas ejected from the 
first generation of stars of GCs can be effectively trapped in the deep gravitational 
potential of their host halos and consequently can be consumed for the formation of the 
second generation of stars without losing a significant amount of gas by ram pressure 
stripping of interstellar and intergalactic medium. During merging of these halos with 
the proto-Galaxy, the halos are completely destroyed owing to the strong tidal field of 
the Galaxy. The self-polluted GCs located initially in the central regions of the halos 
can survive from tidal destruction owing to their compactness and finally become the 
Galactic halo GCs. In this scenario, ejecta of field stars surrounding the central GCs 
can be also converted into stars within their host dwarfs and finally become the second 
generation of stars of GCs. We also discuss the origin of the difference in the degree of 
abundance inhomogeneity between different GCs, such as uj Centauri and NGC 6752, 
in terms of the difference in physical properties between host halos from which GC 
originate. 

Key -words: globular clusters: general - globular clusters:individual (w Centauri)- 
globular clusters:individual (NGC 6752)- galaxies: star clusters - galaxies:evolution - 
galaxies:stellar content 


1 INTRODUCTION 


cesses of stars, such as dredge-up of CN-processed mate¬ 
rial from inner hydrogen-burning regions (e.g., Smith 1987; 
Kraft 1994). Recent observational studies of stellar abun¬ 
dance of some Galactic GCs have revealed star-to-star abun¬ 
dance inhomogeneity among less evolved stars on the main 
sequence, where deep mixing of chemical components are 
highly unlikely (e.g., Cannon et al. 1998 for 47 Tuc). These 
studies accordingly suggested that the primordial pollution 
scenario is more promising than the mixing (or evolutionary) 
scenario (Da Costa et al. 2004; Gratton 2004 for a recent re¬ 
view) . 


The origin of star-to-star abundance inhomogeneity ob¬ 
served in the Galactic globular clusters (GCs) has long been 
discussed based mainly on the following two working hy¬ 
potheses: The primordial hypothesis and the mixing one 
(e.g., Cottrell & Da Costa 1981; Freeman & Norris 1981; 
Smith 1987; Suntzeff 1993; Kraft 1994; Gratton et al. 2004). 
The first hypothesis is that the observed inhomogeneity is 
due to the second generation of stars that were formed from 
gas ejected from the first generation of evolved stars (e.g., 
AGB stars) of GCs (“primordial pollution” scenario, e.g., 
Cottrell & Da Costa 1981). The second is that the observed 
chemical inhomogeneity of GCs can result from internal pro- 


One of key questions related to the primordial pollution 
scenario is whether and how ejecta mainly from AGB stars of 
the first generation of stars can be effectively trapped within 
GCs and consequently converted into the second generation 
of stars. Frank & Gisler (1976) showed that gas of GCs can 
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be efficiently stripped by ram pressure of the Galactic halo 
gas for most GCs. Smith (1996) suggested that the stellar 
ejecta from the first generation of stars are likely to be lost 
entirely from GGs with small binding energies through en¬ 
ergetic outflow of intracluster wind. Gnedin et al. (2002) 
demonstrated that tu Cen could not enrich itself with heavy 
elements of AGB stars owing to efficient ram pressure strip¬ 
ping of the Galactic interstellar medium (ISM), if it formed 
and evolved in isolation. These previous studies thus appear 
to suggest that primordial pollution is not likely within GCs 
evolving in isolation, though numerical attempts have not 
yet been made to investigate the details of primordial pol¬ 
lution processes within GCs. 

Recent numerical simulations have suggested that GCs 
can be formed in the central regions of dwarf galaxies embed¬ 
ded by low-mass dark matter halos (e.g., Bromm & Clarke 
2002; Mashchenko & Sills 2005). Both observational and the¬ 
oretical studies of GCs suggested that massive globular clus¬ 
ters such as uj Cen and G1 were previously stellar nuclei (or 
nuclear star clusters) of nucleated dwarf galaxies (e.g., Free¬ 
man 1993; Bekki & Freeman 2003; Bekki & Chiba 2004). 
Owing to the deeper gravitational potential in the central 
regions of dark matter halos, primordial pollution processes 
can be quite different between GCs evolving in isolation and 
those in the nuclear regions of their host halos. Thus it is 
quite important and timely to discuss whether and how pri¬ 
mordial pollution of GCs can proceed if GCs are within the 
central regions of low-mass dark matter halos that can host 
dwarf galaxy populations. 

The purpose of this Letter is to propose a new scenario 
in which primordial pollution can proceed very efficiently in 
GCs that are located in nuclear regions of their host halos 
at high redshifts. In this scenario, primordial pollution can 
proceed more efficiently in the nuclear GCs than in those 
evolving in isolation, because the ejecta of AGB stars are 
more effectively trapped in GCs (without being significantly 
lost through ram pressure stripping and energetic outflow of 
evolved stars) owing to the deeper gravitational potential of 
their host halos. The nuclear, self-polluted GCs can appear 
as the Galactic halo GCs when the host halos merge with the 
proto-Galaxy and are subsequently destroyed by the strong 
tidal field of the Galaxy (i.e., when field stars and dark mat¬ 
ter halos surrounding the GCs are all removed by the tidal 
stripping). 

We firstly show the three advantages of this scenario in 
explaining why ejecta of the first generation of stars can be 
effectively trapped and converted into stars. Then we discuss 
the origin of the difference in the degree of abundance inho¬ 
mogeneity between GCs in the Galaxy in the context of the 
proposed chemical pollution scenario. Since many authors 
have already discussed advantages and disadvantages of the 
primordial pollution (or enrichment) processes within iso¬ 
lated proto-GC elouds and GCs (e.g., Cayrel 1986; Parmen- 
tier & Gilmore 2001; Thoul et al. 2002; Recchi & Danziger 
2005), we here do not intend to discuss these points. 


2 THREE ADVANTAGES OF CHEMICAL 
POLLUTION WITHIN SUBHALOS 

In the present study, we adopt an AGB pollution scenario 
(e.g., Smith & Norris 1982) in which the second generation of 
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Figure 1. The radial dependences of the escape velocities (Use) 
for four different models: the SB model with Mj, = IO^Mq and 
ro = 0.48 kpc (thick solid), the SB one with ALn = IO^Mq and 
ro = 1.28 kpc (thick dotted), the NFW model with A4h = 10® Mq 
and rs = 0.16 kpc (thin solid), and the NFW one with = 
10®Mq and rs = 0.52 kpc (thin dotted). For comparison, the 
maximum value of the wind velocity of AGB stars is shown by a 
dashed line. 

stars originate from AGB ejecta only after gas ejected from 
Type II supernovae (SNe II) were removed from the systems 
owing to much more energetic processes of SNe II (with the 
wind velocities of an order of ^ 1000 km s“^). Progenitor 
stars of SNe II are considered to have masses of M > lOM© 
(e.g., Woosley & Weaver 1995) and thus have shorter life¬ 
times than AGB progenitor stars with 0.8 < (M/M©) < 8 
(e.g., van den Hoek & Groenewegen 1997). Therefore the 
above assumption on the star formation from AGB ejecta 
only after the completion of SNe II events is regarded as 
reasonable. It is expected that the first and the second gen¬ 
erations of stars are different with each other not in heavier 
elements (e.g., Fe) but in CN-processed ones owing to the 
selective self-enrichment. Although this scenario has several 
problems in explaining guantitatively the observed abun¬ 
dance pattern (e.g., C and O) in GCs (e.g.. Smith & Norris 
1982), we adopt this in order to more clearly demonstrate 
the advantages of the chemical pollution of GCs within 
dwarfs embedded in subhalos. 


2.1 Trapping AGB ejecta 

By investigating whether the terminal (or expansion) veloc¬ 
ity (Uw) of the stellar wind of AGB stars can be smaller or 
larger than the escape velocity (Use) of a halo (or a GC), 
we can discuss whether the ejecta from AGB stars can be 
trapped within the GC for a times scale long enough for 
star formation of the second generation of stars (e.g., Smith 
1996). The observed U, ranges roughly from 2 km s“^ to 40 
km s“^ in Loup et al. (1993) whereas Dupree et al. (1992) 
found evidences of AGB winds with Uv = 90 km s“^. In this 
paper, we discuss whether Use of a dark matter halo can be 
smaller or larger than the possibly maximum value of U, 
= 40 km s“^ (hereafter referred to as Uw,max). In order to 
show the dependences of Use of low-mass halos with differ¬ 
ent masses (Mh) more clearly, we do not intend to derive Uv 
for the combined mass distributions of halos and GCs. 
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For comparison, we adopt the two different profiles of 
dark matter halos at high redshifts: (1) The “SB” profiles 
with large cores proposed by Salucci & Burkert (2000) and 
(2) the “NFW” ones with cnspy cores predicted from the 
standard cold dark matter (CDM) cosmogony (Navarro, 
Frenk, & White 1996). The SB profiles consistent with the 
observed rotation curves of galaxies are described as: 

Psb(r) = , ^ (1) 

(r + ro)(r“ + ro)-^ 

where psb,o and ro are the central dark matter density and 
the core (scale) radius, respectively. For the SB profile, the 
dark matter core parameters, Psb,o, ro, and Mo (where Mo 
is the total dark matter mass within ro) are not free pa¬ 
rameters, and clear correlations are observed between them 
(Salucci & Burkert 2000): 

Mo = 4.3 X 1O^(^)’’^'^M0. (2) 

kpc 

The NFW profile is described as: 

= (r/r.)(l'+r/r.)2’ 

where r, ps, and Ts are the distance from the center of the 
cluster, the characteristic density, and the scale-length of the 
dark halo, respectively. Fig. 1 shows the radial dependences 
of T4sc calculated for the SB model with Mh — lO^M© and 
ro = 0.48 kpc, the SB one with Mh — lO^M© and ro = 1.28 
kpc, the NFW model with Mh — lO^M© and Tb = 0.16 kpc, 
and the NFW one with Mh = lO^M© and rg = 0.52 kpc. 
These values of the NFW models are derived by using the 
lowest mass model in the NFW with a mass-size scaling of 
Mh oc Ta^ reasonable for the CDM halos. 

As shown in Fig. 1, 14sc both for the two NFW models 
are significantly larger than Fw,max in the nuclear regions 
with r/va < 1, which suggests that ejecta of AGB stars can 
be more effectively trapped in nuclear regions of the halos 
compared with their outer ones. The difference between Vesc 
and Vw.max at r/ra < 1 (or r/ro < 1) is larger in the more 
massive model (Mh = IO^Mq) than in the less massive one 
(Mh = IO^Mq) both for the NFW and the SB models. This 
result implies that ejecta of AGB stars can be more effec¬ 
tively trapped in more massive halos and thus that primor¬ 
dial pollution can proceed more efficiently in more massive 
halos. The higher t4sc at a given radius for a given Mh in 
the NFW models in comparison with the SB models sug¬ 
gests that inner density profiles of dark matter halos can 
be an important parameter which determines the details of 
primordial pollution processes. 

Since the mass densities (ph) of halos virialized at red- 
shift 2 : are roughly proportional to (l-|-z)® (Padmanabhan 
1993), t4sc can be larger in halos formed at higher redshifts 
for a given halo mass (i.e., Ksc oc (Mh^ph)^^^ « (1 -|- 2 )^^^ 
for a given Mh). This suggests that ejecta of AGB stars 
can be more effectively trapped in halos formed at higher 
redshifts for a given mass: Primordial pollution of GGs can 
more efficiently proceed in nuclear regions of halos formed 
at higher redshifts. 

2.2 Suppression of ram pressure stripping 

Gas within GGs can be stripped if the following inequality 
is satisfied: 



IMF slope ( —s) 

Figure 2. The dependences of the mass fraction of stars with 
masses ranging from IMq to 8Mq (i.e., progenitors of AGB stars) 
on the slopes of IMF (initial mass function) with —3.5 < —s < 
— 1.5 in GCs for four different models with rrti = O.IMq and 
rriu = 120Mq (thick solid), with mi = 0.5Mq and mu = I 2 OM 0 
(thick dotted), with mi = O.IM 0 and mu = 4OM0 (thin solid), 
and with mi = O.5M0 and mu = 4OM0 (thin dotted). If we 
assume that ejecta of only more massive AGB stars with the 
masses of 6 — 7M0 can be used for star formation, /agb is much 
smaller than those shown in this figure (Bekki & Norris 2005). 


PlSMfrel^ > 2nGT,aT,g, (4) 

where pisM is the density of interstellar gas in halos or disks, 
Wrei is the relative velocity between a GC and the interstellar 
gas, G is the gravitational constant, Eg is the gravitational 
surface mass density, and Eg is the surface density of intr¬ 
acluster gas (e.g., ejecta of AGB stars). By using the for¬ 
mula essentially the same as the above equation (4), Frank 
& Gisler (1976) showed that pisM required for the sweeping 
of gas within the central regions of most GGs is 5 x 10“® 
cm“® for Vrei = 200 km s“^, Es = 10^ M 0 pc“^, and a 
reasonable value of the specific mass loss rate (used for the 
estimation of Eg). The typical mass and tidal radius (rt) 
of GCs are 5 x 1O®M0 and 50 pc, respectively, (Binney & 
Tremaine 1987). This means that mean E^ within rt is by a 
factor of ^ 20 smaller than the above central value and thus 
piSM required for the sweeping of the outer gas is 2.5 x 10“^ 
cm“®. 

It should be noted here that the derived values of pisM 
are much smaller than a typical value of the number den¬ 
sity of H near the Sun (nu ~ lcm“®). The total mass of 
dark matter within the central 50 pc (= n) of the halo in 
the NFW model with Mh = 1O®M0 is 2.8 x 1O®M0, which 
means that Es can be significantly larger in GCs located in 
nuclear regions of the halo than in isolated GCs. This sug¬ 
gests that ram pressure stripping is much less effective in 
sweeping gas of GCs if GCs are located in nuclear regions of 
halos. Although nuclear GCs of halos can keep their gas from 
been stripped by the Galactic halo gas, the surface gravity 
from their host halos appears to be not strong enough to 
prevent gas from being swept by the ISM of the Galactic 
disk with nu ^ lcm“®. 
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2.3 Increase of the total amount of AGB ejecta 


In order to estimate the mass fraction (/agb) of AGB pro¬ 
genitor stars with the masses ranging from IM© to SM© in 
a GC with the total mass of M©, we assume an initial mass 
function (IMF) in number described as ipirni) = 
where mi is the initial mass of each individual star and the 
slope of s = 2.35 corresponds to the Salpeter IMF. The nor¬ 
malization factor A is a function of Md, mi (lower mass 
cut-off), and mu (upper one): 


Mol X (2 - s) 

mu^"® — mi^-® 


( 5 ) 


where mi and mu are regarded as free parameters in the 
present study. 

Fig. 2 shows the dependences of /agb on the IMF slope 
(s) for four different sets of parameters of mi and m^. /agb 
is at most ~ 0.5 in these models. The mass fraction (/ 2 nd) 
of the second generation of stars is described as, 


/ 2 nd oc e* X mej X /agb. 


( 6 ) 


where e* is the star formation efficiency of star-forming gas 
and mej is the mass ratio of gas ejected from AGB stars 
to initial stellar masses of progenitors of AGB stars, e* is 
observed to range from ~ 0.01 (Duerr et al. 1982 for the A 
Sco complex) to ~ 0.4 (Wilking & Lada 1983 for the p Oph 
cloud). If we adopt the table value of mej (=0.46) for AGB 
stars with mi = IM© and Z = 0.001 in van den Hoek & 
Groenewegen (1997), we can find that / 2 nd is at most 0.12 
(= 0.5 X 0.4 X 0.46). 

The derived maximum value of / 2 nd = 0.12 is signifi¬ 
cantly smaller than 0.5 that is required for explaining the 
observed number ratio of CN-strong stars to CN-weak ones 
for NGG 6752 in the context of the primordial pollution sce¬ 
nario (e.g., Smith & Norris 1982). Although / 2 nd required 
for explaining abundance inhomogeneity may well be dif¬ 
ferent between different clusters (thus smaller / 2 nd is still 
acceptable for some GCs), the derived small values of / 2 nd 
imply that if we assume that the second generation of stars 
are formed from the first one initially within GCs, the pri¬ 
mordial pollution scenario can not explain quantitatively the 
observed abundance inhomogeneity of GCs. 

However, if GCs are located in nuclear regions of dark 
matter halos hosting dwarf galaxies, the ejecta of field stel¬ 
lar populations (i.e., major components of the dwarfs) sur¬ 
rounding the nuclear GCs can be fueled into the nuclear re¬ 
gions and converted into new stars there. If a large amount of 
gas ejected from the surrounding stellar populations can be 
converted into the second generation of stars within GCs, 
the mass fraction of the second generation can be signifi¬ 
cantly boosted up: We do not have to assume an extreme 
set of parameter values for IMF and e* in explaining the 
observed / 2 nd (See Smith & Norris 1982 for this problem of 
unusual IMF required for the primordial pollution scenario). 


3 DISCUSSIONS AND CONCLUSIONS 
3.1 Defunct dwarfs as GC hosts 

A key question in this scenario is whether GC host dwarfs 
embedded in dark matter halos can be completely destroyed 
by the strong Galactic tidal field without the central GCs 


being destroyed. Recent numerical simulations have demon¬ 
strated that nuclear star clusters can survive from the strong 
tidal field of the Galaxy owing to their initial compactness 
whereas the main bodies of (nucleated) dwarfs are com¬ 
pletely destroyed and dispersed into the Galactic halo re¬ 
gion (e.g., Bekki & Freeman 2003; Tsuchiya et al. 2003). 
Previous simulations also showed that nucleated dwarfs can 
be more efficiently converted into massive star clusters (i.e., 
naked nuclei) under strong tidal field of groups and clusters 
of galaxies for the SB profiles of dark matter halos (Bekki et 
al. 2003). These theoretical works therefore strongly suggest 
that the present scenario is promising. 

Recent abundance studies of stars both for the Galac¬ 
tic halo and for dwarfs (e.g., dwarf spheroidal, dSph) in 
the Local Group have revealed that [a/Fe] ratios of most 
stars in the dwarfs are generally lower than similar metal- 
licity Galactic halo stars (e.g., Venn et al. 2004). The higher 
[a/Fe] ratios of the Galactic halo can be due to very early 
merging of low-mass dwarf galaxies, which were destroyed to 
form the Galactic halo without efficient chemical enrichment 
of stars (Venn et al. 2004). We propose that these dwarfs 
destroyed in the very early history of the Galaxy (i.e., de¬ 
funct dwarfs) are host galaxies of the Galactic halo GCs with 
abundance inhomogeneity. This proposal is consistent with 
the observed higher [a/Fe] ratios in the Galactic GCs (e.g.. 
Freeman 1993). 


3.2 Diversity in primordial pollution processes 

Smith (1987) divided GCs into three classes according to (1) 
whether star-to-star abundance inhomogeneity can be seen 
for Fe-peak element and (2) whether it can be seen in CN 
abundance (i.e., whether GCs show bimodal CN distribu¬ 
tions). It was found in Smith (1987) that (1) only w Cen 
and M22 show inhomogeneity in Fe-peak elements, (2) GCs 
with no bimodal CN distributions have lower metallicities 
of —2.2 < [Fe/H] < —1.6 (e.g., M92 and M15), and (3) GCs 
with bimodal CN distributions have higher metallicities of 
[Fe/H] > —1.6 (e.g., NGC 6752 and 47 Tuc). However it 
remains unclear why ui Cen shows such inhomogeneity in 
Fe-peak element. 

The present study can provide the following answer for 
the above question, by assuming that the observed inhomo¬ 
geneity is solely due to the primordial pollution processes of 
GCs. uj Cen was formed in the nuclear region of a massive 
subhalo (~ 5 X 10®M©) that was virialized at very high red- 
shift (2 ~ 15) and was more massive than any other subha¬ 
los hosting the Galactic GCs at that redshift. This massive 
halo with a higher mass density (due to earlier virialization) 
could not be soon destroyed by the proto-Galactic tidal field 
because of its stronger self-gravity, and consequently star 
formation could continue for a longer time scale in its nu¬ 
clear region. As a result of this, not only AGB ejecta but 
also some fraction of metal-enriched gas from Type H and 
Type la supernovae could be finally recycled and converted 
into new generations of stars within the u) Gen’s host halo. 
The initially large apocenter distance of its orbit, for which 
a longer time scale of dynamical friction is required for the ui 
Cen’s host to reach the Galactic central region, can be also 
responsible for the prolonged star formation activity (Bekki 
& Freeman 2003). 

Less massive halos containing GCs (e.g., NGC 6752) 
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other than u! Cen had shallower gravitational potential and 
were more easily destroyed by the tidal field of the proto- 
Galaxy in their earlier histories. Only ejecta from more mas¬ 
sive AGB stars with smaller wind velocities (an order of 
~ 10 km s“^) therefore could be converted into new stars to 
become the second generation of stars in GGs without effi¬ 
cient chemical enrichment (i.e., without significant metallic- 
ity spread) due to SNe II with very large wind velocities (an 
order of ~ 1000 km s“^). Less massive halos with smaller 
pericenter and apocenter distances of their orbits would not 
have experienced efficient chemical pollution of GCs owing 
to more rapid destruction of the halos. Thus u) Cen can show 
abundance inhomogeneity both in heavier elements (e.g., Fe) 
and in light ones (e.g., C, N, and O) whereas other GCs can 
show abundance inhomogeneity only in light ones. 


3.3 External pollution scenario 

We have shown that the AGB ejecta can be retained in the 
central regions of their host dwarfs and consequently con¬ 
verted into stars that finally become the second generation 
of stars in GCs. We also have suggested that AGB ejecta 
of field stellar populations surrounding the nuclear GCs can 
be converted into stars to become the second generation of 
stars in GCs. Accordingly, it would be reasonable to say that 
“external pollution” (or modified version of the original pri¬ 
mordial pollution scenario) rather than “self-pollution” is 
a more reasonable jargon that denotes the chemical pollu¬ 
tion processes of GCs in the central regions of dwarfs. This 
external pollution scenario may well explain not only the ob¬ 
served helium overabundance of u) Cen (Bedin et al. 2004) 
but also the large fraction of CN-strong stars in GCs with 
abundance inhomogeneity (Bekki & Norris 2005). 

One of potential problems of this external scenario is 
that the physical mechanism for the selective pollution by 
AGB stars (not by SNe II) is not clearly understood. In or¬ 
der to discuss this point, we plan to investigate gas dynamics 
within the central 0.1 — lOOpc of dwarfs at very high red- 
shifts by using chemodynamical simulations combined with 
the latest results of AGB yields derived by Campbell et al. 
(2005). It is also our future study to investigate whether this 
scenario can explain the observed O-Na and Mg-Al anticor¬ 
relations of GCs. Smith & Norris (1982) suggested that if 
both the first and the second generations of stars are formed 
from ejecta of AGB stars with different mases, the observed 
CN-bimodality in GCs can be explained. This modified pri¬ 
mordial (and external) pollution scenario will be addressed 
by our future chemodynamical simulations in a more quan¬ 
titatively. 

Theoretical studies based on cosmological simulations 
have just started extensive investigation on structural and 
kinematical properties of GC systems of galaxies in a ACDM 
Universe (e.g., Bekki 2005; Yahagi & Bekki 2005). These 
Mpc-scale simulations, combined with sub-pc scale ones on 
star formation within GCs, will provide more robust predic¬ 
tions on spatial distributions and kinematics of the Galac¬ 
tic GCs with different past histories of primordial pollution 
processes and thus be compared with corresponding obser¬ 
vations (e.g., Carretta 2005). 
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